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Abstract . After  a brief  discussion  of  some  undesirable  featixres  of  a number 
of  different  partial  differential  equations  often  employed  in  the  existing 
literature  on  water  waves,  a relatively  simple  restricted  theory  is  constructed 
by  a direct  approach  which  is  particularly  suited  for  applications  to  problems 
of  fluid  sheets.  The  rest  of  the  paper  is  concerned  with  a derivation  of  a 
system  of  nonlinear  differential  equations  (which  may  include  the  effects  of 
gravity  and  surface  tension)  governing  the  two-dimensional  motion  of  incompres- 
sible iiiviscid  fluids  for  propagation  of  fairly  long  waves  in  a nonhomogeneous 
stream  of  water  of  variable  initial  depth,  as  well  as  some  new  results  pertaining 
to  hydraulic  jiamps . 'The  latter  includes  an  additional  class  of  possible  solutions 
not  noted  previously. 
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1 . Introduction 

Kr>,“q at  fitly  it  is  difficult  to  find  t'xnct  solutions  of  cyst«:m£  of  litii  fir  or 
nonlineor  differential  quationa  wiiich  ciifiractorizc  eitiier  a specific  probl'TO 
or  a class  of  problems  in  michanics.  For  tills  reason,  in  many  areas  ol'  mechanics 
of  solids  and  fluids,  widespread  use  has  been  made  in  recent  years  of  approxi- 
mate perturbation  methods,  especially  when  the  systtm  of  differential  '-quatlons 
contains  one  or  more  small  parfimetors . Usually  it  is  not  possible  to  prove 
rigorously  that  the  approximate  solutions  are  related  to  exact  solutions  of  the 
equations  in  a precise  way.  Nevertheless,  considerable  confidence  is  often 
placed  in  the  results,  partly  because  they  may  be  stiown  to  be  satisfactory  in 
special  cases  for  which  exact  results  are  available.  A slightly  different  line 
of  thought  is  to  replace  a system  of  nonlinear  differential  equations  either  by 
an  approximate  linear  system  of  the  same  order,  or  by  a nonlinear  system  of 
lower  order.  In  the  former  case  the  linearization  method  is  clear cut 
and  leads  to  a unique  system  of  equations.  In  the  latter,  the  situation  is 
often  very  confused:  the  procedure  is  singular  in  that  it  reduces  the  order  of 
the  equations  and  may  also  lead  to  some  related  dilemmas  or  questions  concerning 
the  associated  boundary  and  initial  conditions.  Moreover,  the  procedure  is  not 
unique  since  it  may  lead  to  sets  of  apparently  equally  justifiable  equations 
whose  characters  are  quite  different.  This  is  particularly  true  in  the  deriva- 
tion of  nonlinear  differential  equations  which  are  employed  in  the  study 
of  the  propagation  of  fairly  long  water  waves,  as  can  be  seen  from  accounts  of 
the  subject  given  oy  Peregrine  [1]  and  by  Whitham  [2j.  In  fact,  according  to 
V.'hitham  [2,  p.  463],  "The. . .derivations  allow  great  flexibility  and  the 
approach  nat'irally  allows  the  various  alternatives." 

We  list  below  from  [2]  some  of  the  alternative  forms  of  equation.s  for  water 
waves  moving  in  the  direction  of  a fixed  x-axis  for  a stream  of  initial  constant 
depth  h.  Let  the  elevation  of  the  stream  be  h + H and  let  u denote  the  horizontal 


] 


velocity,  where  u,T]  are  functions  of  x and  the  time  t. 
[2,  pp.  1+60-463]  the  system  of  equations 


Furthei',  recall  from 


Tl^+  {u(h-f  11)]^  = 0 , 

u^tuu^fgTl^.i  c^hTl^^x  = ° 


and  the  pair  of  equations  attributed  to  Boussinesq,  i.e., 


(1) 


T]^+  [u(htT])]^  = 0 


u.  + uu  + 

t X 


^ 3 "\tt 


= 0 


(2) 


wh<re  subscripts  indicate  partial  differentiation  with  respect  to  t or  x , 

2 

c = gh  and  g is  the  acceleration  due  to  gravity.  Both  systems  of  equations  (l) 
and  (2)  allow  for  wave  propagation  in  either  direction  along  the  x-axis.  For 
waves  moving  along  the  positive  x-direction  only  there  is  the  Korteweg-deVries 
(hereafter  referred  to  as  the  K.dV.)  equation  [3]  given  by 

° (3) 


or  the  following  equation  due  to  Benjamin  et  al.  [4]: 

Some  of  the  properties  of  the  differential  equations  (l)  to  (4)  should  be 
noted.  It  may  inmediately  be  verified  that  the  set  of  equations  (2)  and  (4) 
have  steady  state  solutions  only  if  T]  and  u are  both  constants.  Also,  although 
the  K.dV.  equation  (3)  admits  a solitary  wave  in  which  the  velocity  at  infinity 
is  zero  when  the  stream  there  is  at  its  undisturbed  height  h,  it  does  not  admit 
a steady  state  solution  with  u constant  and  1^  = 0 at  infini+y.  This  fact  is 
related  to  another  property  of  (3)  which  is  also  shared  by  (2)  and  (4):  the 
three  sets  of  equations  (2)  to  (4)  are  not  invariant  in  form  under  a cotistant 
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superpored  rigid  body  moiion  of  the  whole  fJuid.  To  see  this,  suppose  that  e 
constant  supei'posed  rigid  body  trnns  1 at.ions i velocity  k is  imposed  on  the  wliole 

*f  f 

fluid  so  that  the  particles  at  the  place  x ai-e  lisplaced  to  x at  time  t where 

X = X + kt  , t = L t- k , (k  constant)  , ('^) 


and  the  velocity  u is  replaced  by  u + k.  'Iiien,if  the  dependent  variables  x,t  in 
(2)  to  (4)  ore  changed  to  (5),  the  orpiations  for  u,Tl  in  terms  of  x^,t^  are  dif- 
ferent from  those  in  terms  of  x,t.  For  example,  equation  (3)  becomes 


T1  +k'n^  + c(n||)T]^-Hjch%^^  +=0 


(6) 


t X x XXX 

Tills  means  that  the  character  of  the  solutions  of  (2),  (3)  and  (4)  are  radically 
altered  by  superposing  a constant  rigid  body  translational  velocity  on  the  fluid, 
which  is  contrary  to  what  happens  if  we  use  the  full  tiu’ee-dimensional  equations 
of  motion  for  an  inviscid  fluid.  On  the  other  hand,  the  set  of  equations  (l) 
are  not  subject  to  this  drawback,  and  they  do  have  useful  steady  state  solutions. 
It  may  be  argued  that  because  of  the  nature  of  the  approximation  in  obtaining 
(2)  to  (4)  from  the  three-dimensional  theory  we  should  not  expect  these  equations 
to  be  invariant  under  a superposed  constant  translational  velocity,  but  this  then 
leaves  in  doubt  which  version  of  any  of  the  sets  (2)  to  (4)  are  correct  and  are 
to  be  chosen  as  basic.  'Ihe  difficulty  disappears  if  we  linearize  any  of  the 
above  sets  since  the  resulting  equations  are  then  invariant  under  a small  super- 
posed constant  translational  velocity,  as  we  would  expect. 

It  might  appear  from  the  above  discussion  that  the  equations  (l)  may  be 
preferable  to  any  of  (2)  to  (4),  but  arguments  are  put  forward  in  [2,  p.  462] 
to  suggest  that  the  set  (2)  is  to  be  preferred  to  (l).  Although  considerable 
use  has  been  made  of  some  of  the  equations  (1)  to  (4),  it  would  seem  that  they 
all  rest  on  a somewhat  shaky  physical  foundation. 

In  recent  years  the  derivation  of  nonlinear  equations  which  are  suitable 
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for  the  propagation  of  fairly  long  water  waves  has  been  approached  from  a 
completely  different  point  of  view,  tDimcly  via  the  tlieory  of  a (iirected  (or 
Cosserat)  surface.  A direct  two-dimensional  theory  of  this  kind  is  established 
in  [5j6]  for  water  waves  employing  Integral  balance  laws  for  mass  conservation, 
momentum,  moment  of  momer.tiun  and  (when  required)  energy,  in  parallel  with  what 
is  done  in  the  development  of  the  three-dimensional  theory,  and  is  valid  for 
water  of  variable  initial  depth.  Exich  equations  are  fully  nonlinear  and  satisfy 
all  usual  invariance  requirements  as  do  their  three-dimensional  counterparts. 

The  resulting  field  equations  of  this  theory  contain  certain  assigned  fields 
and  inertia  coefficients  which  are  unspecified,  but  these  are  identified  by  an 
appeal  to  rather  easily  accessible  results  which  can  be  deduced  from  the  tlu"ee- 
dimensional  equations.  An  alternative  derivation  given  in  [7]  starts  from  the 
tliree-dimensional  energy  equation  and,  after  employing  a single  approximation 
for  the  velocity  field,  proceeds  only  with  the  help  of  invariance  conditions 
under  superposed  rigid  body  motions  and  leads  to  the  same  system  of  equations 
as  those  in  [5 >6] . 

The  present  paper  is  concerned  with  some  extensions  of  previous  work  [5,6] 
by  a direct  approach.  First,  we  briefly  construct  a restricted  theory  of 
directed  sizrfaces  in  which  some  of  the  kinematic  and  kinetic  variables  are 
restricted  at  the  outset.  'This  restricted  theory,  which  may  be  regarded  as  a 
special  case  of  a more  general  theory  of  a Cosserat  surface,  is  particularly 
useful  in  application  to  problems  of  fluid  sheets  and  is  utilized  in  the  rest 
of  the  paper  for  incompressible  fluids  in  which  the  mass  density  may  vary  with 
depth.  A nonlinear  system  of  (two-dimensional)  differential  equations  is 
derived  for  propagation  of  water  waves  in  a nonhomogeneous  stream  of  variable 
initial  depth  [see  h)qs.  (35)>  (5^)  and  (55)]  and  some  extensions  of  earlier 
work  [6,  section  8]  on  hydraulic  jumps  are  discussed. 
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G«"nfral  lifickj/i-oimd . A Ro»;trictrd  TT-icory  of  Directed  r.ijri'accs. 


We  first  discuss,  in  ttiis  fi'ction,  some  bnckj^round  information  from  the 
tiieory  of  a Cosscrat  (or  directed)  siirface  and  then  siunmarize  a special  case  of 
the  theory  in  a mnnner  wiiicii  is  particularly  suitable  for  application  to  the 
propagation  of  water  waves  in  a noniiomuge-neous  fluid.  We  recall  that  a Cosserat 
surface  C comprises  a material  surface  (embedded  in  a Eiaclidean  3-space)  and  a 
single  deformable  director  attached  to  every  material  point  (or  particle)  of  C- 
Let  the  particles  of  the  material  surface  of  C be  identified  with  a system  of 
convected  coordinates  9 (a  =1,2)  and  let  the  S’j'face  occupied  by  the  material 
surface  of  C in  the  present  configuration  at  time  t be  referred  to  by  J.  Let 
r and  d denote  the  position  vector  of  a typical  point  of  J and  the  director  at 
the  same  point,  respectively,  and  also  designate  the  base  vectors  along  the 
8°^-curves  on  by  a . Then,  a motion  of  the  Cosserat  surface  is  defined  by 
vector-valued  functions  which  assign  position  r and  director  d to  each  particle 
of  C at  each  instant  of  time,  i.e.f 


r = r(e“,t)  , d = d(e“,t)  , [a, ad]  > 0 


(7) 


ind  the  condition  (7)-,  ensures  that  the  director  d is  nowhere  tangent  to  J. 

3 ~ 


The  base  vectors  a and  their  reciprocals  a , the  unit  normal  a^  and  the 

~Qf  ~ 

con^jonents  of  the  metric  tensors  a and  a®^  at  each  point  of  J are  defined  by 

op 


dr 

^ ‘ a? 

i 

.2c 


„0  -O 

> a . a-  = , a 


= a 


crp  ^ 


QfP  Of  3 
= a • « 


(8) 


® ~3  Si  ^ Se 


, a = det  a^  , a^  = [a^^^a^]  > 0 , 


^Recall  that  when  the  particles  of  a continuum  are  referred  to  a convected 
coordinate  system,  the  numerical  values  of  the  coordinates  associated  with 
each  material  point  remain  the  same  for  all  time. 

^Ihe  choice  of  positive  sign  in  (7)3  is  for  definiteness.  Alternatively,  it 
will  suffice  to  assume  that  [a^^dj  ^ 0 with  the  luiderstanding  that  in  any  given 
motion  the  scalar  triple  product  Toj^d]  is  either  >0  or<0. 
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where  6^  is 
are  given  by 


the  Kronecker  delta. 


The  velocity  and  the  director  velocity  vectors 


V = ^ , w = d , (9) 

where  a superposed  dot  denotes  differentiation  with  respect  to  t holding  0® 
fixed.  Throughout  this  paper,  we  use  standard  vector  and  tensor  notations. 

In  particular,  Greek  indices  take  the  values  1,2  and  the  usual  summation  con- 
vention over  a subscript  and  a superscript  is  employed. 

A general  theory  of  a Cosserat  surface  by  Green,  Naghdi  and  Wainwright  [8] 
is  developed  within  the  framework  of  thermodynamics;  the  derivation  in  [8]  is 
carried  out  mainly  from  an  appropriate  (two-dimensional)  energy  equation, 
together  with  invariance  requirements  under  superposed  rigid  body  motions. 

Here  we  adopt  the  mode  of  derivation  of  the  basic  theory  given  by  Naghdi  [9] 
or  the  more  generalized  fom  of  the  theory  employed  by  Green  and  Naghdi  [6]. 
Special  cases  of  the  general  theory  can  be  obtained  by  introduction  of  suitable 
constraints,  thereby  resulting  in  constrained  theories.  Alternatively,  cor- 
responding special  cases  can  be  developed  in  which  the  kinematic  and  the  kinetic 
variables  are  suitably  restricted  a priori  and  then  restricted  theories  are 
constructed  by  direct  approach.  Such  special  cases  of  the  general  theory  have 
been  discussed  previously  by  Naghdi  [9>  Secs.  10  and  I5]  and  by  Green  and 
Naghdi  [10]  and  are  of  particular  interest  in  the  context  of  elastic  shell 
theory.  In  the  present  paper,  however,  we  construct  anotlier  type  of  restricted 
th(.-ory  which  is  particularly  suitcKi  for  application  to  problems  of  flxiid  sheets. 
'Ibe  resulting  equations  can  also  be  obtained  as  a constrained  case  of  those  given 
for  directed  fluid  sheets  in  [6],  but  we  find  it  more  convenient  to  restrict  the 
kinematic  and  the  kinetic  variables  at  the  outset  and  construct  a corresponding 
restricted  theory  from  an  appropriate  set  of  conservation  laws  in  integral  form. 
We  therefore  confine  attention  here  to  a theory  in  which  the  director  d,  while 
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d'''formlng  nlong  its  length,  always  remains  parallel  to  a fixed  direction 
specified  by  a constant  tinit  vector  b.  Tims,  I’ocalling  al.so  (il),  we  have 

d = 0(9“, t)b  , w = w(e“,t)b  , w = 0 . (10) 

Let  P,  bounded  by  a closed  ciorve  9P,  be  a part  of  J occupied  by  an  arbitrary 
imiterial  region  of  the  surface  of  c in  the  present  config\iration  at  time  t and 
let 


V = v°^a  = V (11) 

~ ~a  CK~ 

be  ' tinit  normal  to  &P.  It  is  convenient  at  tide  point  to  define 

cei  aditional  quantities  as  follows:  The  mass  density  p=p(0^,t)  of  the 

sixrface  J in  the  present  configuration;  the  contact  force  N=N(0'^,t  ; v)  and 
the  contact  director  force^  M=M(0^,t  ;v),  each  per  unit  length  of  a curve  in 
the  present  configuration;  the  assigned  force  f=f(0^,t)  and  the  assigned 
director  force  £=£(9^,t),  each  per  unit  mass  of  the  surface  V;  the  intrinsic 
(surface)  director  force  m per  unit  area  of  J',  the  inertia  coefficient  k=k(0''^) 
which  is  independent  of  time  and  is  associated  with  the  director  velocity;  the 
specific  internal  energy  e = e(0^,t);  the  heat  flux  h = h(0"^,t  ; v)  per  \init  time 
and  per  unit  length  of  a curve  5P;  the  specific  heat  supply  r = r(0^,t)  per  unit 
time;  and  the  element  of  area  do,  and  the  line  element  ds  of  the  siirface  J. 
F’irther,  in  view  of  the  assumed  form  (lO)  for  the  director,  we  express  M,  m 
and  A in  terms  of  their  components  along  and  perpendicular  to  the  unit  vector 
b,  i-e., 


^The  terminology  of  director  couple  is  also  used  for  M depending  on  the  physical 
dimension  assumed  for  the  director  d.  Here  we  choose”  d to  have  the  physical 
dimension  of  length  so  that  M has  the  same  physical  dimension  as  N.  For 
further  discussion  see  [6].  ~ 
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M = M(9^,t  ; v)b  + b X S(9'^,t  ; v)  , S-b  = 0 , 

m = m(9'^,t)b-tbxs(9^,t)  , s*b  = 0 , (12) 

^ + b X c(9'^,t)  , c • b = 0 , 

where  M,  m and  I are  scalar  functions  and  S,  c are  vector  fimctions  of  their 

ar^’:uments . Also,  it  is  convenient  to  decompose  the  assigned  fields  f and  i 
into  two  parts,  one  of  which  represents  the  three-dimensional  body  force  acting 
on  the  continuum  which  are  assumed  to  be  derivable  from  a potential  function 
D(r,ii)  and  the  other  which  represents  the  effect  of  applied  surface  loads  on 
the  major  surfaces  of  the  fluid  sheet.  Thus,  we  write 

£ = - f ♦ £ ’ « = (-  M * b . (13) 


With  the  foregoing  definitions  of  the  various  field  quantities  and  with 
reference  to  the  present  configuration,  the  conservation  laws  for  a restricted 
theory  of  a Cosserat  surface  (different  from  that  discussed  in  [9]  or  [10])  are: 


_d^ 

dt 


dt 


r p do  = 0 , 

V 

;^|pvda=|Pfdo+|  Nds  , 

~ ~ J^p~ 

b^j*  pkwda  = b[[  (pji-m)do  + f MdE]+bx[|  ( pc-s  )do  + j*  S ds  ] , 

~dtjp  ~ Jp  ~ Jp 

[ prxvdo=  [ p[rxf  + dx(bxc)]do+[  [rxN-tdx(bxS)]ds  , 
at  ,]p~  ~ .)p  ~ ~ ~ ~ ~ ~ ~ ~ ~ 

f p[ e + Ci  + |-( V • V + kw^ ) ]do  = | p(r  + f • v + £w )do  + j (N  • v 4 Mw  - h)ds 
Jp  ~~  Jp  ~~ 


(li+a) 

(I4b) 

(l4c) 

(I4d) 


(I4e) 


In  the  above  equations  (l4a)  is  a statement  of  conservation  of  mass,  (l4b)  the 
conservation  of  linear  momentum,  (l4c)  that  of  the  conservation  of  the  director 
momentum,  (l4d)  the  conservation  of  moment  of  momentum,  and  (l4e)  represents 
the  conservation  of  energy.  It  should  be  noted  that  the  quantities  M and  £ 
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make  no  conti-itmt ions  to  the  moment  of  mumentum  equation,  and  the  quantities 
o and  r make  no  eontri but j on  to  tlie  energy  equation  in  tlie  prcs'-nt  lu-.stri ote-d 
theory. 

■'nde»’  suitable  continuity  a.sstjmptions,  tiie  ciirve  force  N,  the  director 
force  M iind  the  heat  flux  h can  be  expressed  as 


N = N V , 

~ ~ a 


M = m“v  , S = s“v  , S“  • b = 0 , 


h = qv^  , 

where  the  fields  are  functions  of  0'^,t.  The  five  conservation  equa- 

tions in  (l4)  then  yield  the  local  equations 


pa^  = y(e'^)  , 


(a^N^)  + yf  = yv  , 

~ jO  ~ ~ 

1 jy  1.1  1 

(a^Nr)  ^ + YJi  = + Ykw  , (a^s“)  ^+yc  = a^s  , 

a xN“+dx(bxs)  + d x(bx  s“)  = 0 , 

pr  - q*^i  - pc  + . V + mw  + = 0 , 

|o  ~ ~5Ck  ,a 

where  a comma  denotes  partial  differentiation  with  respect  to  the  surface 
coordinates  0^  and  a vertical  line  stands  for  covariant  differentiation  with 
respect  to  the  metric  tensor  of  the  surface  J.  It  should  be  noted  that  the 
vector  fields  S®  and  s are  workless  and  do  not  contribute  to  the  energy  equa- 
tion (20). 
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3 • Incompressible  Inviscid  Fluids. 


P 


Kithin  the  scope  of  the  ttieui-y  of  n Cosscrot  surf'ucc,  pi’i  viously  in  [6] 
constitutive  e-quations  have  been  obt.nineci  for  a fluid  sheet  wiiich  model  the 
properties  of  the  three-dimensional  isothermal  inviscia  fluid.  r"or  the 
present  restricted  theory,  we  first  note  that  the  conditions  of  inctjmpressibility 
reduce  to  the  single  condition^ 


^ = 0 . 

With  the  help  of  (9)  and  (10),  the  condition  (2l)  can  be  rewritten  as 


(21) 


.op  , 

l>e  a.  X b 


V wa_  • b 

~3  ~ 


0 , 


(22) 


where  e ^ is  the  alternating  tensor.  For  an  incompressible  inviscid  fluid  at 
constant  temperature,  the  response  functions  N®,m,M“  are  workless,  i.e.. 


• V + mw  + M^w  = 0 

~ ~,q;  ,a 


(23) 


provided  v ,w  satisfy  the  constraint  condition  (22).  It  can  then  be  shown  that^ 

~,Q' 

n“  = - , m = - p^a^  • b , = 0 , (24) 

where  p^  is  an  arbitrary  scalar  function  of  e'^jt.  Moreo.er,  in  a manner  similar 
to  that  indicated  in  the  appendix  of  [6],  for  an  inviscid  fluid  at  constant 
temperature  we  also  have 

q“  = 0 , e = 0 . (2b) 

It  then  follows  from  (20)  that 

* 

I ” - ■ , ■ ■ . ■ , 

) d 

, In  general,  there  are  two  conditions  of  incompressibility  in  the  theory  of 

• incompressible  directed  fluid  sheets  given  in  [6).  In  our  present  development, 

{ since  d is  specified  by  (lO)^,  the  second  condition  is  satisfied  identically  and 

1 I the  corresponding  pressure  (arising  lYom  the  constraint  response)  is  a part  of 

' the  response  functions  for  C®  and  s. 


^The  development  leading  to  (24)  is  similar  to  a more  general  discussion  in 
section  4 of  f6l. 


r = 0 . 


(26) 


Also,  for  the  present  restricted  tiieory  the  moment  of  momentum  equation  (19) 
is  satisfied  if 


I s0  +•  p 0c*^^a„(h  • a ) = 0 


(27) 


Tile  last  equation  is  consistent  with  the  restrictions  s“  . b = 0,  s • b=0  given 
in  (15}  and  (12).  From  (18)  and  (27),  after  eliminating  s,  we  obtain  the 
following  equation  for  S : 


1-  .1  Q 

(0a^S‘^)  +p  0a^e®^a. (b  • a ) + y0c  = 0 

^ 5 Of  O <^.<P  ^ rs^ 


(28) 
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4 . Water  Waves  for  NonhomoKcneous  Stream  of  Variable  Depth. 

Our  previous  discussion  of  water  waves  in  [6]  was  confined  to  onc-diniensioi.al 
flows.  Here  we  consider  two-dijnensional  flows  and  at  the  same  time  allow  our 
model  to  reflect  the  (three-dimensional)  properties  of  a nonhomogeneous  incom- 
pressible fluid.  Let  e^  ,e„,e_  be  a set  of  right-handed  constant  orthonormal 

~1  ~3 

base  vectors  associated  with  rectangular  Cartesian  axes  and  let  the  position 
vector  r in  (7)^^  and  the  director  d in  (lO)^  be  represented  as 

r = xe  + ye  + ijte  , d = 0 e , b = e , (29 ) 

~ ~ ~ ~3 

1 2 

where  x,y,i)i,0  are  functions  of  9,0  ,t.  The  velocity  v and  the  director  velocity 
now  take  the  forms 

V = ue,  + ve.  + \e  , w = we_  , (30) 

r*<i 

where 

u=x  , v=y  , , W=0  (31) 


and  we  note  that  the  velocity  components  u,v,\,w  may  be  regarded  as  functions 
1 2 

of  either  9 ,0  ,t  or  of  x,y,t.  From  (30)  follow  the  expressions 


and 


ue,  + ve_  + \e_ 


(32) 


u = + UU  + VU  , V = V,  + uv  -(  vv  , 

t X y ’ t X y ’ 


\=\^  + uK^  + v\y  , w = w^  + uw^+vWy  , 


(33) 


where  as  in  section  1 the  subscripts  x,y,t  designate  partial  differentiation 
with  respect  to  x,y,t  when  u,v,X,w  are  regarded  as  functions  of  x,y,t.  Also, 
by  (7),  (8)  and  (29),  the  base  vectors  a and  the  unit  normal  a can  be  expressed 
in  the  forniB 


12. 


+ e t 


(iX  . ^ 1 . 

a - — r e + —V  ^ ^ f , 

ae^  ae^'  ae^  ''^ 


ax 


ae 


2 XI 


ae 


. ay  aj 

«■,  + — ^ e„  + — ^ e , 

ae“ 

-ALu^I  ■ 

a „2s 


~3  a(0^9^ 

With  the  use  of  (30)  and  (33)  > tlie  incomprecsibility  condition  (22)  assiames 
the  simpler  form 


(34) 


S(u  +v  )+w  = 0 
^ X y' 

and,  with  the  help  of  (24),  (32)  and  (33),  the  equations  of  motion  (I7)  and 


(18)^ reduce  to  the  forms 


yu  = yf  ' e 


a(e\e^)  ^ 


yv  = y£  • ^ 


. a(x,y) 


1 .2,  ^y  > 


a(e  ,0  ) 


yx  = y£  • 


ykw  = yl  + 


a(x,y). . E 
a(9\9^) 


where  p is  given  by 


(35) 


(36) 

(37) 

(38) 
09) 


P = P^0 


(40) 


It  rtmiainE  to  specify  values  for  the  coefficients  Y,k,  the  assigned  force  f 
and  the  assigned  director  force  1.  For  this  purpose,  we  are  guided  by  the  cor- 

responding  fluid  sheet  in  the  t^lree-dimensional  theory  in  which  an  incompressible 

♦ * 

homogeneous  fluid  under  gravity  -g  e^  fJows  over  a bed  specified  by  the 
position  vector 


*We  use  g*  (instead  of  g)  for  gravity  and  reserve  the  symbol  g for  later  use 
in  Bq.  (48). 
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and  wo  specify  the  surface  of  the  fluid  by 


p = + y^  I- B(x,y,t)e^  . (42) 

In  (4l),  O'  is  a giv>’n  Ainction  of  x,y  but  3 in  (42)  depends  on  x,y,t.  At  the 
surface  (42)  of  the  stream  there  is  a constant  pressure  p^  and  a constant 
normal  surface  tension  T.  At  the  bed  the  (unknown)  pressure  p depends  on  x,y 
and  t.  Thus,  the  fluid  moves  with  the  surface  (42)  and  the  normal  pressure  p 
at  this  surface  is 

P*'  = Po-'l  ’ (’^3) 


where 


q = P — pITp • 

(l+p2  2 3/2 

At  the  bed  (4l)  the  normal  velocity  of  the  fluid  is  zero  and  the  pressure  p 
takes  the  value 


(44) 


* 

P 


p(x,y,t) 


(45) 


where  p is  to  be  determined. 

To  proceed  further,  we  introduce  a set  of  convected  coordinates  0^  (i= 1,2,3), 

3 

let  the  surface  0=0  coincide  with  the  surface  and  consider  the  three- 
dimensional  region  of  space  between  the  surfaces  (4l)  and  (42)  occupied  by  the 
fluid.  Any  point  in  this  three-dimensional  region  is  then  specified  by 

We  also  suppose  ’liat  the  siurfaces  0^  = co  - g and  0^  = u)+2>  with  m a constant  to 
be  detennined  later,  are  coincident  with  (4l)  and  (42),  respectively,  so  that 


I 


a - 1(1  I (ou  - 2 '>0  < ^ - iji  • (oi  ^ 2 )■/) 


The  nasF  J'.'iiRiLy  p of  the  fluid  ic  amun'^d  to  vary  witli  d'-'pth  uo  that 

p p ; ) • 

'ITie  base  vectoi-t-  in  tin'  fibm'-e-mfnt i oned  tlir  i ’-d ijU' nsi cria  1 r'Uii'n.  y: 
w 11  ;is  their  conjugated  and  j-elat-d  "'in  h'.-  calculated  fr  .m  (hu) 

b\it  we  only  j'ecord  hei'e  the  fotinulna 


1 P ’ 

a(e  ! 


= '-U'  ' - ('V  ' 9'^^s.,)eo  " 

X X y y Me  ,B  ) 


V,(  now  make  use  of  the  results  derived  in  [6]  in  order  to  obtain  explicit  values 
of  y,f,£  and  the  parameter  <«.  First,  in  relation  to  the  top  and  bottom  surfaces 

3 

of  the  fluid,  we  choose  the  surTaca  9=0  so  that  the  center  of  mass  of  the  (three- 
dimensional)  fluid  region  under  consideration  a j ways  lies  on  this  s'urface  and  we 
then  identify  this  surface  with  the  siarface  ^ in  the  tlicoi-y  of  Coss-Tat  surface. 
This  leads  us  to  impose  the  condition 

_y  I 

^p'(0^)0^de^  = 0 (-Huj<2)  , (49) 

(JO  - ;>- 

which  determines  y,>.  Tlien,  recalling  the  results  in  [6],  we  have^ 

4.  1 4-  1 

pa‘"  = V = J ^P  g^de^  , = J ^(B"')  p g^d9  (90) 


or  equivalently 


where 


V = ^0  , Yk  = K„0  , 

3(9  ,3  ) 3(0  ,9“^) 

f ^ 4-1 

30-’  , <2  " J ^ • 


^fee  particularly  equations  (3-7)  and  (3-19)  'f*  [■■>]. 


15. 


Also,  witl'  the  use  of  (13)>  (5l)  find  (48),  from  Fqs.  (3-15)  of  [6]  for  the 
assigned  force  and  the  assigned  director  force  we?  obtain 


v£  = U(Po-q)0x  - ^ 

(q-P  ^p-g  <0)e^J  ■ , 

° d(e\e^) 

yl  = [(q-P  )(|  + ou)  i-p(-2-  + uu)]  ■ , 

° d(9\9^) 

YC  - [ (Pj^-q)(s  +uj)3j^' p(-y  + 

' ^(e  ,9  ) 

- r(p„-q)(2-  + u))3„-p(-|  ^'iu)Oyl£p 
ri  = g i|f  • 


(53) 


Substitution  of  (51)  to  (53)  into  (36)  to  (39)  results  in  the  differential 
equations  of  mo + ion 


K0U=-p^+  (p^-q)0^-po^  , 

K0v=-Py+  (p^-q)By-Pay  , 

• — * 

K7»X  = q - Pq  + P “ g > 

= (q-Po)(|  + 'Ju) +P(-2 +uj) +1  • 


(54) 


f-loreover,  since  the  bed  of  the  stream  is  stationary,  from  (4?)  and  (SB)^ 
have 


Of  = uof  * va  = i {til  - 1)0  = \ *■  (iu  - ^)v  . (55) 

A y 

The  aoove  system  of  equations  is  independent  of  the  rfjnaining  equations  (2?) 
and  (28)  which  involve  S®,s.  The  fields  S®,s  correspond  to  appropriate  con- 
straint responses  for  the  restricted  motion  (lO)  considerpd  in  the  pr^-sent 
paper; and,  since  they  are  not  completely  determined  by  (27)  and  (28),  there  is 
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^ some  arbitrariness  in  ttieir  specification.  We  do  not  need  to  consider  ei^uations 

(27)  and  (28)  further  lure. 

• To  E\imm£irize  the  results  obtained  so  far,  we  observe  that  the  propagation 

I of  water  waves  in  a nonhomogeneous  Incompressible  inviscid  fluid  is  governed  by 

the  system  of  differential  equations  (35)  j (5^0  and  (55)  with  the  coefficients 
’''■><2  given  by  (52)  and  (u  determined  from  (49).  Wlien  the  fluid  is  liomogencous 

4 

with  p a constant,  it  readily  follows  from  (49)  and  (52)  that 

uu  = 0 , iC=p  , *^2"12  ■ 

With  values  (56),  the  differential  equations  (35) 5 (54)  and  (55)  become  identical 
with  those  derived  previously  [7]  from  the  tliree-dimensional  theory  of  homogeneous 
incompressible  inviscid  fluids.  For  unidirectional  flow  along  the  x-direction, 
the  equations  reduce  to  those  obtained  in  [6]  by  a direct  approach. 

The  above  equations  have  been  derived  from  the  integral  balance  equations 
(l4a-e).  With  the  help  of  (24),  (25),  (26),  (30),  (5l)  and  (53),  these  integral 
balance  Laws  can  be  reduced  to  the  forms; 


-j^K0dxdy  = O , 

u dxdy  = f [ (p^-q)e  - pa  Idxdy  - [ p dy  , 
at  Op  jp  o X X j^p 

I KJv  dxdy  = f [ (p^-q)  P,  - por  Idxdy  + \ p dx  , 
at  jp  Jp  o y y j^p 

j*  K ^ A.  dxdy  = r (q-P.+P-g  K^)dxdy  , 
at  ..p  Jp  o 

^ f dxdy  = + (q-p^^)(2-  + 'ju,!  - p('2  -'u)]dxdy  , 


‘i . Solitary  Waves  on  n Stream  with  Level  Bed. 

Wh'.-n  tVie  bed  of  the  stream  Is  li.;vel  and  the  fluid  is  homogeneous  with  a 
constant  value  for  p , it  iias  been  shown  in  f5]  that  the  system  of  r-quatjons 
(3!?)j  (5^)  and  (55)  admits  a solution  in  the  form  of  a solitary  wave.  A similar 
solxition  is  possible  here  even  when  the  IMuid  is  nonhomogeneous  with  the  mass 
d-  nsity  p varying  with  depth.  (Imitting  details,  for  a wave  travelling  along 
the  x-axis  with  a speed  c,  the  solitary  wave  has  the  form 


— = 1 + p,  sech  [ ( P — 

(l-2(a)<t4K, 


c^p,(x-ct) 
ch 


-] 


(63) 


where 


M-  = 


(c-  - <)-a- 


gh(l-2(ij)  5 c >1 


(64) 


Previously  [5]  the  notation  \ was  employed  for  a quantity  corresponding  to  p, 
defined  by  (64)2-  the  present  paper,  the  symbol  X is  utilized  for  a 
different  purpose  in  (35)- 


6.  Hydraulic  Jiunps. 


In  a previous  paper  [6,  Pec.  8]  wo  discussiil  hydraulic  Jumps  at  a step  in 
the  horizontal  bed  of  a stream.  Here  we  extend  the  discussion  to  nonhomofjeneous 
fluids  in  which  the  mass  density  varies  with  depth,  and  at  the  same  time  indicate 
the  nature  of  an  additional  class  of  possible  solutions  wliich  was  not  noted  in  [6]. 

We  assume  steady  flow  of  the  fluid  parallel  to  the  x-direction  in  a stream 
whose  bottom  is  level  except  for  a finite  Jump  at  the  origin,  and  thus  we  specify 
or  by 


(u)  - = a = 


0 for  X < 0 , 

d for  X > 0 , 


'li„  = (5  - uu)0^  for  X < 0 , 


(65) 


where  d is  a constant.  We  neglect  surface  tension  so  that  q = 0,  for  convenience 
set  the  surface  pressure  p^=0,  and  consider  only  steady  motions  on  either  side 
of  X = 0 with  the  possibility  that  the  wave  height  0 at  x = 0 has  a stationary 
finite  discontinuity.  At  x = 0,  we  assume  that  0,u,X.,w,p  change  from  the  values 
0^,u^,X.^,w^,p^  to  the  values  utilize  the  one- 

dimensional form  of  the  equations  (57),  (58),  (60),  (61)  and  (62)  with  v = 0 in 
a usual  way  to  obtain  the  appropriate  jump  conditions.  Thus^ 


U202  = = n 5 

(66) 

Kn(u2-u^)  = Pj^  - P2  +X  , 

(67) 

= X2  = (2-0u)w^  = (?-U))w2  , 

(68) 

K*  -M- 

-Uj)g  nK(0^-0^)+g  nKd  = ‘ P2^^2  ' 

(69) 

^The  constant  n in  (66)  to  (69)  corresponds  to  k in  the  corresponding  jump 
conditions  given  in  [6,  Sec.  8]. 
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In  [6,  S>^c.  8]  the  ciiscussion  was  limited  to  the  special  case  In  which  X = 0, 
However,  in  gi-nernl , the  quantity  X could  be  nonzero;  it  arises  from  the  inte- 
grated value  par  in  (^8)  since  p or  may  become  large  at  the  step,  p being  the 
pressure  on  tiie  bottom  of  the  stream.  The  term  X on  the  right-hand  side  of 
(t;/),  tli'T’-fore,  represents  ti*e  resultant  force  exerted  by  the  step  on  the 
fluid  measured  along  the  positive  x-directionf  Assuming  that  tiie  fluid  rfjmains 
in  contact  with  the  step  tliroughout  the  motion,  we  impose  the  conditions 


X S 0 when  d > 0 , 

X § 0 when  d < 0 


(70) 


In  order  to  illustrate  the  differences  which  arise  when  X^O,  we  consider 
now  the  simplest  problem  discussed  in  [6,  Sec.  8].  Thus,  when  the  stream  is 
level  on  either  side  of  the  step,  from  (5^)  we  have 


"^  2 n 

01  = h » = (i -U))g  Kh  , ^ , 


<t>2  ~ ^ > ?2  = (?-u))g  ) ^2  H ’ 


(71) 


For  convenience,  we  introduce  the  notations 


T ^ n 

L — , , T — ^ ^ 

^ g (l-2,u)h3 


(72) 


and  then  substitute  (71 ) into  (67)  and  (69)  and  also  use  (66)  to  obtain 


(73) 


Kg  h 


’^Similar  resultant  forces  occur  in  jump  conditions  utilized  by  Caulk  [11]  in 
his  treatment  of  the  problem  of  fluid  flow  under  a sluice  gate  based  on  the 
two-dimensional  theory  of  directed  fluid  sheets  [5»6]. 


A 
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To  Gijnplify  the  dificussion  tlitit  follows,  v;e  set 


4t^  = T ^ , ftTg  - - T ♦ (t^+8t)^  , 


(74) 


T,  +i-  = (4-  + 2t)'^ 


and  ttien  rewrite  (73 )i  o 

I 


± ^ (.  -l)(l^t  ) , 


(7^) 


- ^ (1-2uj)(t  -L)(L+t  +1)  . 

Kg  h " 

If  we  confine  attention  to  values  of  t>1,  then  Tj^>I  and  1 for  all  values 
of  lu- 

In  the  rest  of  this  discussion,  we  confine  attention  to  homogeneous  fluids 
in  which 


(D  = 0 , K = p , To+5  = (-4  + 2t)2  (t>1,  t >1)  , 


(76) 


but  a parallel  discussion  can  be  carried  out  for  nonhomogeneous  fluids  with 
(ju/O.  The  previous  solution  [6,  Sec.  8]  corresponds  to  the  case  for  which 


L = 


d 


= ^ ° ’ h = "irr 


(77) 


On  the  other  hand,  if  we  take  the  value  of  d to  be  that  specified  by  (77)^} 
then  in  addition  to  (77)^  we  find  a second  possible  solution  for  L,  namely 

T^  + 1 T^+1  5 T.,+1  1 

which  satisfies  (75)^^.  It  can  be  verified  from  (7^)2  second  value  of 

L given  by  (78)  corresponds  to 


X > 0 , d > 0 


(79) 


in  comformlty  with  (70). 
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To  contin\ie,  we  observe  that  for  values  of  d/h  specified  by 


0<i< 


(t^-1  )' 


(80) 


there  are  in  general  two  positive  values  of  L which  satisfy  (75)^  but  only  one 
of  these  is  compatible  with  (79) • Purther , if 


:t3-d- 


- h h 


(81) 


where  d^  is  a maximum  positive  value  of  d for  varying  L in  (75)^>  then  there 
are  two  values  of  L corresponding  to  each  value  of  d/h  and  both  of  these 
satisfy  (79)- 

A more  general  discussion  of  the  other  problems  discussed  in  [6,  Sec.  8] 
can  be  given  in  a similar  manner,  but  we  do  not  pursue  the  matter  further. 
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